Note on the Origin of the Highest Energy Cosmic Rays. 
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In this note we argue that the galactic model chosen by E.-J. Ahn, G. Medina- Tanco, P.L. Bierman and 
T. Stanev in their paper discussing the origin of the highest energy cosmic rays, is alone responsible for the 
focusing of positive particles towards the North galactic pole. We discuss the validity of this model, in particular 
in terms of field reversals and radial extensions. We conclude that with such a model one cannot retreive any 
directional information from the observed direction of the cosmic rays. In particular one cannot identify point 
sources at least up to energies of about 200 EeVfJ] . Therefore the apparent clustering of the back-traced highest 
energy cosmic rays observed to date cannot be interpreted as an evidence for a point source nor for the identification 
of M87, which happens to be close to the North pole, as being such a source. 



It was recently claimed W that the Ultra High 
Energy Cosmic Rays (UHECR) observed up to 
now, may originate from a common source lo- 
cated in the direction of the North pole of our 
Galaxy; the M87 cluster (Virgo) being the best 
candidate for such an emitting object. This con- 
clusion was drawn from the result of the back- 
tracing of the 13 highest energy cosmic rays ob- 
served to date (assumed to be protons, or helium 
nuclei for the most energetic ones) in a modeled 
galactic magnetic field. In their study the field 
originates from a "galactic wind" (analogous to 
the solar wind) extending far away from the visi- 
ble part of the galactic disk. 

In this paper we study the mathematical prop- 
erties of such a field, and show that the observed 
convergence is an intrinsic property of the field 
model, and not an evidence of a unique, pointlike 
source for the highest energy cosmic rays. 

1. Focussing properties of the proposed 
magnetic wind 

As far as ultra-high energies (over 100 EeV) are 
concerned, the bending effects are dominated by 
the long range behaviour of the galactic field. In 

1 EeV for Exa electron Volts or 10 18 eV 



the model proposed in the asymptotic field, 
in spherical coordinates, is purely azimuthal and 
reads : 



where the normalization constant BqTq = 
70 /LtG.kpc is defined from the local value of the 
field in the Solar system. The bounds of the re- 
gion where the galactic wind extends are not well 
defined; the authors use 1.5 Mpc in their numer- 
ical simulations. 

The most important feature of such a field (in 
the absence of a cutoff on r) is that the bending 
integral J B v dr is divergent in any radial direc- 
tion except the polar axis. As a result, whatever 
the energy, a charged particle can never escape to 
infinity in a direction other than a pole. In prac- 
tice, this strong focussing effect is limited by the 
cutoff, which therefore plays a crucial role. Using 
the field and radial limits given above the bending 
integral is about 500 x sin 6 EeV along a radial 
trajectory. In other words particles of 100 EeV 
will only escape within a cone of less than about 
10 degrees around the polar axis. 

When considering particles of a given sign the 
orientation of the radial component of the Lorentz 
force depends on the polar projection of the ve- 
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Figure 1. Trajectory of protons of 100 to 160 EeV, back-traced from the Earth's northern (left) and 
southern (right) hemispheres, the Galactic North pole is to the right. Top : Transverse coordinate 
projected on the galactic meridian plane containing the Earth versus range. Bottom : Polar angle versus 
range. 



locity, therefore the "positive" pole (as defined by 
the orientation of curl B) is focussing, while the 
"negative" pole is antifocussing. 

These features are intrinsic to the field model 
(especially its slow decrease with the distance) 
and we can suspect that the convergence of the 
trajectories found in Jl[] is therefore not related 
to any specific property of the observed set of 
highest energy cosmic rays. 

2. Numerical simulations with random 
data 

To confirm this hypothesis, we have numeri- 
cally back-traced, in the field model of Ref. |l]], 
a random set of cosmic rays drawn from a uni- 
form distibution on the (Earth) sky. As most 
of the observered high energy events are actually 
concentrated around 10 20 eV, and since the au- 
thors of [|J have further assumed that the two 
highest energy events could be Helium nuclei the 
range of magnetic rigidity of the observed data 
sample is quite narrow. For a fair comparison 
we have therefore generated a sample of protons 
with a flat energy distribution ranging from 100 
to 160 EeV. 



As expected the trajectories drawn on Fig.[j] 
clearly show the strong focussing effect of the 
field. Most of the focussing take place over the 
first few 100 kpc as mentioned in (!]]. We have 
drawn separately the trajectory for particle reach- 
ing the Earth's northern hemisphere and reaching 
the southern one. The focussing effect is stronger 
and faster for rays originating from the nortbQ 
which is expected given the configuration of the 
Earth's rotation axis with respect to the galactic 
center. Using a random population we obtain the 
same behavior as in |lj showing that the focussing 
is a property of the field and not of the events. 

3. Discussion and conclusion 

First the validity of the model may be ques- 
tioned: if there are reversals of the azimuthal 
component of B in the galactic disk (as acknowl- 
edged in 0), how can the wind remain consistent 
with a "coherent" parametrization sin6/r at long 
distances (several 100 kpc, much more than the 
distance between regions with reversed fields) ? 

2 The 13 events used in |l| were all observed in the north- 
ern hemisphere, which is where the largest detectors are 
installed. 
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Figure 2. A representation similar to the plot of 
Ref Black (open) circles are for rays hitting 
the Earth's northern (southern) hemisphere. 



Normally one would expect some destructive in- 
terferences between the wind contributions com- 
ing from diffcrents parts of the disk, hence an in- 
tensity decreasing faster than 1/r; the argument 
that most cosmic rays are observed in the direc- 
tion opposite to the galactic center (where no re- 
versal occurs) is not valid, because their bending 
depends mainly on the long range behaviour of 
the field. 

If however the model of Ref.JIJ] is true, then 
the accumulation of events at the pole is not at 
all an evidence for a pointlike source of the ob- 
served rays. This model only demonstrates that 
our sensitivity on extragalactic charged particles 
might be limited to a small solid angle (decreasing 
with energy) around the galactic polar direction, 
whatever their initial origin could be. 

One should note that even if the extragalac- 
tic flux is isotropic the integrated luminosity at 
Earth would be the same with or without this 
galactic wind. Despite the strong dispersion in 
the original directions the restriction of the angu- 
lar acceptance due to the focussing effect is com- 
pensated by the collecting area. 



Addressing the question of the active galaxy 
M87 (Virgo A) as a possible source of UHECR, 
the only possible affirmative conclusion is: if the 
field model of [Q is valid, and if the sources are 
known pointlike objects, a possible candidate is 
M87. However, as acknowledged in jjj, this sce- 
nario implies a regular transverse magnetic field 
from here to Virgo, i.e. over a length of about 
20 Mpc, of about 2 nanogauss. Therefore the 
overall system would behave like a spectrum an- 
alyzer : a magnetic spectrometer (the transverse 
field) followed by a collimator (the Galactic field), 
strongly selecting the initial momentum of the 
cosmic rays. 
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